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In the present paper, we investigate the kaon twist-3 distribution amplitudes (DAs) cf>^ a within 
the QCD background field approach. The SU / (3)-breaking effects are studied in detail under a sys- 
tematical way, especially the sum rules for the moments of (f>p a are obtained by keeping all the mass 
terms in the s-quark propagator consistently. After adding all the uncertainties in quadrature, the 
first two Gegenbauler moments of are a^ p (lGeV) = -0.376t^, a^ p (lGeV) = 0.70lig:^, 
a^ CT (lGeV) = -0.160±^4 and a| f <T (lGeV) = 0.369t^g, respectively. Their normalization 
parameters /i^|iGcv = 1. 188±g'g|§ GeV and /x^liGcv = 1.021+^055 GeV - A detailed discussion on 
the properties of p>cr moments shows that the higher-order s-quark mass terms can indeed provide 
sizable contributions. Furthermore, based on the newly obtained moments, a model for the kaon 
twist-3 wavefunction ^!^ a (x, kx) with a better end-point behavior is constructed, which shall be 
useful for perturbative QCD calculations. As a byproduct, we make a discussion on the properties 
of the pion twist-3 DAs. 
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I. INTRODUCTION 

Meson distribution amplitude (DA), which describes 
the momentum fraction distribution of the parton in me- 
son, is an important component for the QCD light-cone 
sum rule (LCSR) and the QCD factorization theory [J- 
5]. In dealing with the exclusive processes, it is conve- 
nient to arrange the meson's DA by its different twist 
structures. The leading-twist DA shows the momentum 
distribution of the valence quarks in the meson, which 
usually provides major contribution to the QCD exclu- 
sive processes. The higher-twist DAs describe either the 
contributions from the higher Fock states with additional 
gluons and / or quark-antiquark pairs or the contribu- 
tions from the transverse motion of quarks (antiquarks) 
in the leading-twist components. Usually, the contribu- 
tions from the higher-twist DAs are power suppressed to 
that of the leading-twist in the large Q 2 -region. How- 
ever, the twist-3 DAs may provide sizable contributions 
for certain cases, so it arouses people's more and more 
interests, c.f. Refs. 

Kaon twist-3 DAs are important input parameters for 
the kaon electromagnetic form factor, the B — > K tran- 
sition form factor and etc., whose properties have been 
investigated within the QCD sum rules and the kx fac- 
torization approach accor dingly [IM1]. More precise 
data are coming at LHC, it would be useful to study 
the higher-order and higher-power suppressed contribu- 
tions so as to provide a deeper understanding of standard 
model parameters. For example, it has been pointed out 
that the SU f (3)-breaking effect is about 10% for B -> K 
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transition form factors, so a careful study on K meson 
distributions shall lead to a better estimation of these 
form factors. 

We have studied the QCD sum rules for the pionic 
twist-3 DAs in Ref.[ll[, which are based on the frame- 
work of the QCD background field theory [H|-[24| . In the 
present paper, we shall improve our technology adopted 
there and then investigate the twist-3 DAs of K meson 
by carefully dealing with its S f J7/(3)-breaking effect. Such 
an effect is responsible for the different behaviors between 
kaon and pion DAs. 

Basic assumption of the QCD sum rules is the in- 
troducing of the nonvanishing vacuum condensates such 
as the quark condensate (qq) and the gluon condensate 
(G 2 )- Different to the conventional SVZ sum rules [2o| . 
the background field approach provides a systematic de- 
scription for these vacuum condensates from the view- 
point of field theory. And, it is convenient to derive use- 
ful relations among different non-perturbative matrix el- 
ements. Under the background field approach, it assumes 
that the quark and gluon fields are composed of the back- 
ground fields and the quantum fluctuations around them. 
Nonperturbative effects can be described by the vacuum 
expectation values of these background fields, while the 
calculable perturbative effects are expressed by quantum 
fluctuations. Then to take the background field theory as 
the theoretical foundation for the QCD sum rules, it not 
only has distinct physical picture, but also can greatly 
simplify the calculation due to its capability of adopting 
different gauge conditions for quantum fluctuations and 
background fields respectively. 

Because of the influence from background fields, the 
quark and gluon propagators shall include nonperturba- 
tive component inevitably. For the SVZ sum rules, one 
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usually takes the following quark propagator formula [2t 



S(x, 0) 



i $ 
2ir 2 x' 



12 



(m) 



192 



(g s qaTGq) 



1152 



(g s qaTGq) $ 



4ir 2 x 2 
■■, (1) 



where • • • stands for even higher-dimensional terms and 
higher-order mass terms. Note that the above quark 
propagators in configuration space are given as an ex- 
pansion in quark mass, and the mass terms are kept only 
up to hrst order. For the light-quark propagators, it is 
enough. However, the omitted higher-order mass terms 
may lead to sizable contributions to the meson or baryon 
with heavy quark(s). Even for the case of K meson, the 
contributions from higher-order s-quark mass terms, ei- 
ther positive or negative, are sizable. Hence to obtain 
a better understanding of S'C//(3)-breaking effect for K 
meson, one needs to take the sizable mass terms into 
consideration in a more proper way. 

The remaining parts of the paper are organized as fol- 
lows. In Sec. II, we present the calculation technology for 
deriving the sum rules for the moments of the kaon twist- 
3 DAs. And a model for the kaon twist-3 wave functions 
is also presented. Numerical results are given in Sec. Ill, 
where the properties of kaon twist-3 DAs are discussed. 
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Sec. IV is reserved for a summary. In the Appendix, we 
give useful formulas for simplifying the matrix elements 
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II. CALCULATION TECHNOLOGY 
A. Sum rules for the pseudoscalar twist-3 DAs 

Under the background field theory, the quark and 
gluon propagators satisfy the following equations [H| : 



(i p -m)S F (x,0) =5 4 (x) 
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- ^ -iffTM" and (-D M ) = 6 ab d p - gf abc A<^ 
are gauge covariant derivatives in the fundamental and 
adjoint representations respectively, 
ground gluon field, G° v = d^A^ - d v A« , .,, .,, 
is the gluon field strength tensor and f abc (a,b,c = 
1, 2, • • • , 8) is the structure constant of SU (3) group. Fix- 
ing the gauge freedom of the background field by the 
fixed-pointed gauge [13, HI], x^A® = 0, we obtain 
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for the quark propagator and 
d 4 a 
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for the gluon propagator, where b^ vp = |G !//1 (0) and 



= §[G^ ;p (0)+G w;i ,(0)]. Here G v „{x) 
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g s T a G° (x), the gauge invariant function G VKp \ 

g s T a Dp b G b ^(x)\ x= o, and the symbol ••• stands for the 
irrelevant terms for our present analysis that will lead to 
higher-order operators over dimension-six. To use a prop- 
agator in momentum-space form as Eq. Q has been sug- 
gested in the literature already, e.g. it has been suggested 
to deal with the D*Dn and B*Bn couplings in Ref . [29jj . 
However in these discussions, usually the first two terms 
in the quark-propagator are kept only. For the present 
case, one may find that the third term should be kept to 
provide a more accurate sum rules up to dimension-six 
operators. As a special case, by taking only the first- 
order mass term, we can obtain the quark propagator in 
the coordinate space, 
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Since the quark propagator in momentum space keeps the 
mass terms naturally, so we shall adopt (j4j other than ([6]) 
to do the following calculation. In fact, as will be shown 
later, the high-order mass terms are indeed important for 
giving a more sound 5?7/(3)-breaking effect in 4>p a . 
The pseudoscalar twist-3 DAs <f>^ and 4>a are defined 
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where f = 2x — 1, qi = d and q 2 — u for pion, q\ = s 
and q 2 = u for kaon, respectively; the parameters fp 
and [iPp" stand for the decay constant and the normal- 
ization parameter of the pseudoscalar, respectively. The 
DA moments are defined as 

{Q) p = [ dx(2x - (7) 
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which satisfy 
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respectively. In deriving the sum rules for the moments, 
we adopt the following correlation functions: 
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Fig-© shows the Feynman diagrams for the pseu- 
doscalar DA moments, where the background gluon fields 
are included in the Fermion propagators implicitly and 
the background quark fields are depicted as crosses. 



FIG. 1. Feynman diagrams for the pseudoscalar DA moments, 
where the background gluon fields are included in the Fermion 
propagators implicitly and the background quark fields are 
depicted as crosses. The left big dot stands for the vertex 
75(12 • D n ) and cr MI/ 75(12 • Z)(" +1 )) for (frp and respectively. 



Following the standard QCD sum rule technology, the 
sum rules of the moments can be derived. And we obtain 
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Because the current quark mass of u(d)-quark is quite 
small, we have set m 2 = m u ~ 0, mi = ~ for 
pion and mi = m s for kaon accordingly. M stands for 
the Borel parameter, mp is the pseudoscalar mass, and 
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trix elements: (a s G 2 ) = (a s G A v G A ^), (g s qaTGq) = 
(g s qa^T A G A ^ q ), (g 3 JG 3 ) = (glf^G^G 3 ""^) 
Q. If setting n = in Eqs. (|llll2|l . one can obtain the 
sum rules of the normalization parameters. In deriving 
the sum rules (|11I12[) . we have implicitly adopted the 
following Borel transformation formulas 
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and we have used the simplified matrix elements 
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detailed derivations are presented in the Appendix. 

From the sum rules (|11|12[) . it is found that their 
perturbative parts and the dimension-four gluon con- 
densate part come from Fig.(QJi); the dimension-three 
quark-antiquark condensate part and the dimension- 
five quark- gluon condensate part come from Fig.([T|3); 
and the dimension-six four-quark condensate part comes 
from Figs.(QJ>nJn}i). Numerically, it can be found that 
the contribution for (f)^ from Fig.([TJ;) is very small (~ 
m 2 /M 6 ). But if one adopts the quark propagator ((6|), 
the contribution of (f>p from Fig. (JTJ;) should be twice than 
that of Fig-dlji) (~ 1/Af 4 ) and is sizable. This shows 
that by keeping the mass-terms properly, one can obtain 
a correct estimation of the relative importance among 
different Feynman diagrams. To show the s-quark mass 
effects more clearly, we shall discuss the different conse- 
quences caused by the using of the propagators (|4]) and 
© in the next section. And we shall find the importance 
of using the propagator Q , which keeps the higher order 
mass-terms in a more consistent way. 

It is well-known that the kaon twist-3 DAs can be ex- 
panded in Gegenbauler polynomials as 



CO 

^( M/ , x) = 1 + J2 a n K J^ f )C 1 J 2 (2x - 1), 



(13) 



1 Note in the condensate (g s qq) 2 , whose the coupling constant g s 
comes from the gluonic background field, so we should treat the 
condensate as a whole. 



a K P (^f) an d a Ka(^f) are Gegenbauler moments at the 
factorization scale (if. 

With the help of the above sum rules (0 [8|) for the DA 
moments, we can obtain the kaon twist-3 DAs, e.g. 



a K,p 



l K,i 



3 (£p)#r' a K,p - y (i€) K 3) 



K 



l K,cr 



2 

35 
12 



1 



K 



(15) 



And the Gegenbauler moments at any scale (i can be 
obtained from the renormalization group equations from 
an initial factorization scale [16f : 
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where Tn = C F (l - ( » +1) 2 (n+2) + 4 Yj£l 
(lliV c - 2JV/)/3 with C F = 4/3. 
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B. kaon twist-3 wavefunctions 

The kaon wave function and its DA can be related with 
the following equation, 



<l>p,A x i Vf) = 
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where (if ~ 0(lGeV) is the factorization scale. Due to 
the renormalization group equation (|16l) , the distribution 
amplitudes under different choice of (if ~ 0(lGeV) can 
be related with each other through evolution, which shall 
result in the same behaviors at the present considered 
accuracy [3(j. Hereafter for definiteness, we set (if = 
lGeV. 

Following the same idea of Refs.QS EE M, 
where its transverse momentum dependence is con- 
structed on the BHL-prescription [34[, the kaon twist-3 
wave functions can be constructed as 
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TABLE I. The normalization parameters and the first two moments of 4>p,<j- 



M 2 (GeV 2 ) 
04 ) 2 (GeV 2 ) 


[0.875, f.f03] 
1.287- f.5f7 


A/ 2 (GeV 2 ) 


[0.800, 1.198] 
-0.126 ±0.010 


A/ 2 (GeV 2 ) 

m« 


[1.088, 1.094] 
0.425 ± 0.001 


M 2 (GeV 2 ) 
(^^)(GeV 2 ) 


[0.968, 1.034] 
1.178 - 1.247 


M 2 (GeV 2 ) 


[1.103, 1.267] 
-0.094 ±0.001 


M 2 (GeV 2 ) 


[0.825, 1.117] 
0.329 ±0.013 
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•,3/2, 



exp 
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m 9jS indicate the constituent quark mass, and their stan- 
dard values are m„ ~ 0.30 GeV and fh s ~ 0.45 GeV. 
The parameters A^ a , Bp a , C^ a and (3^ a can be deter- 
mined by the average value of the transverse momentum 

(ki>^ (« (0.350GeV) 2 (H), 
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the wave function normalization 
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and the first two DA moments 
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K meson masses = 139.57018 ± 0.00035 MeV and 
rriK = 493.677 ± 0.016 MeV; the pion and kaon decay 
constants / w = 130.41 ± 0.20 MeV and f K = 156.1 ± 0.8 
MeV. The vacuum condensates have been calculated and 
updated since 1979 EL c.f. Refs.[H, We take 

the dimension-four and dimension-six condensates to be 
0: (a s G 2 ) = (7.5 ± 2.0) x 10~ 2 GeV 4 , ( 5s 3 /G 3 ) = 
(8.3 ± 1.0)GeV 2 x (a s G 2 ). And for the quark condensate 
and quark-gluon condensate we take [4l[ : (uu) (2GeV) = 
-(0.254±0.015) 3 GeV 3 , (dd) = (uu), (ss) / (uu) = 0.74± 



0.03, (g s qaTGq) 



qq) with ml = 0.80 ± 0.02GeV^ 



and (g s qqf = (2-7±g;|) x lO^GeV 6 . The continuum 
threshold parameter is taken to be around the mass 
square of the first exciting state of the meson. Consider- 
ing the first exciting states are 7r(1300) and iT(1460) for 
pion and kaon respectively [36|, we take s P T ' a — 1.69±0.10 
GeV 2 and s P jf = 2.13 ± 0.10 GeV 2 . The leading order 
a s is fixed by a s {M z ) = 0.1184 ± 0.0007 and the 
renormalization scale is taken as M. 



B. 
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III. NUMERICAL ANALYSIS 
A. input parameters 

From the Particle Data Group [36|, we take the cur- 
rent s-quark mass as m s (2GeV) = lOOj^o MeV; 7r and 



To derive proper Borel windows for the sum rules of 
<$£ a (x, /if), the criteria are to suppress the unwanted 
continuum contribution and the higher-dimensional con- 
tribution as much as possible so as to obtain more accu- 
rate results. 



First, we determine the normalization parameters jA p £ . 
In Refs. 0, EH , it is calculated by using the idea of the 
quark equation of motion (QEM). While it has been 
pointed out that the quarks inside the meson is not ex- 
actly on-shell [13], so the results in Refs. 0, HI} is only an 



approximation. As a notation, it is found that the three- 
particle twist-3 distributions ^-k^k can be related with 
(t>p a through the QEM H3 . However due to the simi- 
lar reason, we do not discuss the three-particle distribu- 
tions with those relations in the present paper. A simple 
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M 2 (GeV 2 ) 



FIG. 2. The first and second moments of kaon twist-3 DAs versus the Borel parameter M 2 , where the shaded bands are the 
uncertainties caused by varying all the input parameters within their reasonable regions. 



TABLE II. Main uncertainties for the normalization parameters and the first two moments of <fip l<T , where the uncertainty of a 
particular parameter is obtained by fixing other parameters to be their center values. 





Oi- ) 2 










&)« 




-0.014 
+0.003 


-0.077 
+0.044 


-0.045 
+0.029 


-0.025 
+0.014 


-0.038 
+0.024 


-0.009 
+0.003 


gPf 


+0.049 
-0.051 


+0.050 
-0.053 


±0.000 


±0.014 


±0.000 


+0.007 
-0.008 


(a 3 G 2 ) 


±0.041 


±0.039 


T0.002 


±0.044 


to.ooi 


±0.033 


(ss) 


-0.004 
+0.003 


-0.023 
+0.019 


-0.011 
+0.009 


-0.013 
+0.011 


-0.017 
+0.015 


-0.020 
+0.017 


(uu) 


+0.018 
-0.016 


=F0.001 


-0.011 
+0.010 


+0.012 
-0.011 


±0.000 


-0.001 
+0.000 


(g s soTGs) 


±0.000 


+0.002 
-0.001 


±0.001 


+0.006 
-0.005 


+0.008 
-0.006 


+0.022 
-0.018 


p.CT 






+0.008 
-0.010 


-0.027 
+0.033 


+0.006 
-0.009 


-0.021 
+0.032 










-0.017 
+0.014 





discussion on this point can be found in Ref. 14|, where 
compatible results for <g; a 3 ^> and / 37r with those de- 
rived from QEM [6] have been obtained through proper 
consideration. At the present, by setting n = in the 
sum rules (|11|12[) . we can obtain the sum rules for fj^f ' . 
To set the Borel window for pPjf , we take the continuum 
contribution to be less than 40%, and the dimension- 
six condensate contribution to be less than 2% for \F K 
and 4% for [i a K . The values of {p? K ) 2 and (m^tMk) an ^ 
the their corresponding Borel windows are collected in 
TablH which is obtained by setting all the input parame- 
ters to be their center values. Main uncertainties caused 
by the current quark mass m s , the continuum threshold 
s P jf , the dimensional operators and etc. are collected 



in Tab|Tl] Other smaller uncertainties caused by the pa- 
rameters as (g s uaTGu), (g s ss) 2 , (g s uu) 2 , fx and etc. 
are not presented. By taking all uncertainty sources into 
consideration, we obtain, 

M^IigcV - l-188i°;°43 GeV, (24) 
MkIigcV - 1.02lt°;°j» GeV, (25) 

where the renormalization group equation of \jP k and fi K 
[44l [45[ has been adopted to run its value from the scale 
M to 1 GeV. Note our value of /^IigcV is different 
from the value obtained by the on-shell condition (i.e. 
A^rliGoV - 1-424 GeV [H}) by about (17 ±3)%. 



Second, we calculate the first two moments of kaon 
twist-3 DAs. The Borel window for the second moment 
of (£,p a ) K is determined by setting the continuum con- 
tribution to be less than 35% and the dimension-six con- 
densate contribution to be less than 3%. The Borel win- 
dow for the first moment of (£,p a ) K is determined by 
setting the continuum contribution to be less than 5% 
for (£a) K and to be less than 3% for {^) K ; while the 
dimension-six condensate contribution is set to be less 
than 1% for and less than 3% for (&) K . The re- 

sults together with the corresponding uncertainties are 



presented in TabU and TabHIl To show the uncertain- 
ties more clearly, we draw the first two moments of the 
kaon twist-3 DAs versus M 2 in Fig@, where the shaded 
bands are the uncertainties caused by varying all the 
input parameters within their reasonable regions. By 
adding these uncertainties in quadrature, and with the 
help of the relation between different moments (TT51) and 
the scale running relation (|16l) . we can obtain the cor- 
responding Gegenbauler moments at the scale /z/ = 1 
GeV: 

< p (lGeV) = -0.376±S, (26) 
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TABLE III. Wavcfunction parameters for ipp (x, kx) with typical DA moments at /if = 1 GeV. 



a K,p 
2 

a K,p 


-0.524 
0.210 0.701 1.182 


-0.376 
0.210 0.701 1.182 


-0.273 
0.210 0.701 1.182 


Bp 
o p 
/3p A '(GeV) 


163.047 139.472 126.284 
-0.249 -0.188 -0.137 
1.133 1.559 1.968 
0.506 0.549 0.581 


168.526 143.459 129.680 
-0.063 -0.019 0.021 
1.168 1.583 1.986 
0.503 0.546 0.578 


172.755 146.721 132.586 
0.064 0.096 0.129 
1.194 1.601 1.999 
0.501 0.544 0.575 


TABLE IV. Wavefunction parameters for ip^ (x, kj_) with typical DA moments at /if — 1 GeV. 


2 

a K,a 


-0.234 
0.220 0.369 0.532 


-0.160 
0.220 0.369 0.532 


-0.109 
0.220 0.369 0.532 


Af(GeV" 2 ) 

c? 

ft? (GeV) 


194.004 150.888 116.789 
-0.080 -0.062 -0.049 
0.114 0.225 0.349 
0.442 0.479 0.519 


199.642 154.708 119.355 
-0.007 0.010 0.026 
0.129 0.237 0.360 
0.440 0.477 0.516 


203.549 157.461 121.467 
0.043 0.059 0.075 
0.138 0.246 0.367 
0.438 0.475 0.514 



4 jP (lGeV) 
a* (lGeV) 



0.701 



+0.481 
-0.4911 



-0.160 



+0.051 
-0.074' 



(27) 
(28) 



a 2 Ka (lGeV) = 0.369 



-0.163 
-0.149- 



(29) 



Based on the above moments and the formulas pre- 
sented in Sec. II, we can obtain the kaon twist-3 wave 
function parameters, which are collected in Tabs. (|IIIHV[) . 
Here Tabs. (IIIIHV[) correspond to the factorization scale 
(if = 1 GeV. With the help of Eq. ((T7|) , we can obtain 
the kaon twist-3 DAs <fi^ and cj>^ , which are presented 
in Fig.©. Our <p^ (Left diagram) or (Right dia- 
gram) are drawn by the solid lines which are defined by 
Eq.(fP7"|). and the dash-dot lines are for Eqs. (|13ll4[) with 
Hf = 1 GeV, respectively. As a comparison, we also give 
the DAs of Ref . [13| under /if — 1 GeV, which are drawn 

I 



I 

by the dashed lines. Here, in doing the comparison, we 
need to replace ^^(x) in Ref. [3] to be <^ CT (1 — x), be- 
cause in Ref.[i~3j x stands for the momentum fraction of 
s-quark; while in the present paper, x is taken as that 
of u (or d) quark. These two figures indicate that the 
kaon twist-3 DAs, especially cj>p , have a better end-point 
behavior. Such a BHL-improved behavior shall be help- 
ful to obtain a reasonable result for kaon related pro- 
cesses, such as the kaon electromagnetic form factor and 
the B — >• K transition form factor with kx factorization 
approach or LCSR. Some previous calculations can be 
found in Refs.llTlIlElal. 



TABLE V. Contributions of the higher order quark mass terms to various vacuum condensate parts for the moments of 4> PtlT . 
For estimation of the contribution of higher order quark mass terms to every vacuum condensate, the percentage of which is 
obtained by calculating the ratio of higher order mass terms (m n , n > 2) before each vacuum condensate with those of (m™, 
n < 1). 













(a s G 2 ) 


4.3% ~ 2.8% 


2.3% ~ 2.2% 






m q (qq) 


-0.7% ~ -0.4% 


17.6% - 17.5% 


1.2% ~ 1.0% 


2.8% ~ 1.9% 


m q (g 3 qaTGq) 


7.2% ~ 4.4% 


-8.0% 7.9% 


11.7% ~ 10.0% 


-0.6% ~ -0.4% 


(g s qq) 2 


-16.5% 9.1% 


-0.5% 


-1.6% ~ -1.3% 


1.0% ~ 0.7% 
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FIG. 3. Kaon twist-3 DAs <j>g(x,Hf) (Left) and <j>*(x,(if) 
(Right), where the solid and the dash-dot lines are for our 
DAs defined by Eq.fTf)! and Eqs. (|13ll4|) with ju = 1 GeV 
respectively. As a comparison, the DA of Ref.( 13j) at fif — 1 
GeV is shown by a dashed line. 



As has been argued in the Introduction in order to 
provide a more sound estimation on the S'L r /(3)-breaking 
effect in the if-meson involved processes, we need to use 
the much more complex Eq.(j4|) other than Eq.([6]) as the 
quark propagator. To show this point clearly, we show 
in Tab ]V] how the higher order mass terms contribute to 
the corresponding vacuum condensates for the first two 
moments of <f>p a . For estimation of the contribution of 
higher order quark mass terms to every vacuum conden- 
sate, e.g. (a s G 2 ), m q (qq), m q (g s qaTGq) or (g s qq) 2 , the 
percentage of which is obtained by calculating the ratio of 
higher order mass terms (m™, n > 2) before each vacuum 
condensate with those of (m n , n < 1). Tab |V] indicates 
that because the s-quark mass is not small, it shall lead 
to sizable contributions. For example, its contribution to 
m q (qq) for K can be up to 17%. 

IV. SUMMARY 

The background field approach provides a systematic 
description for the vacuum condensates from the view- 
point of field theory and it provides a convenient way 
to derive the QCD sum rules. We have made an in- 
vestigation over the kaon twist-3 DAs a within this 
approach. Furthermore, the SUf (3)-breaking effects are 
studied in detail under a more systematical way, espe- 



cially the quark propagator (Q} that keeps the mass terms 
consistently is adopted. As have been shown by TabfVl 
higher-order mass terms can indeed provide sizable con- 
tributions to the kaon DA moments. For example, its 
contribution to m q (qq) for (%p) K can be up to 17%. So 
to obtain a sound estimation for the S'C//(3)-breaking ef- 
fect, we need to take these higher-order mass terms into 
consideration. Moreover, such a propagator shall also be 
helpful for deriving information on the meson or baryon 
with heavy quarks. Some more works on its application 
to the heavy meson/baryon properties are in progress. 

As for the kaon twist-3 DAs <f>p and <p^ , we have 
studied their normalization parameters and moments 
within the QCD sum rules under the background field 
approach. For its normalization parameters, we obtain 
MkIigoV = GeV and ^ K \ 1GcV = 1.021±g;g|g 

GeV. As for the moments of 4>p a , around /// — lGeV, 
we obtain ($) K = -0.126l Q ;°| , K = 0.425™, 

(a) K = -0.094±°;°3° and (&) K = 0.329l°; O og. Basing 
on these moments, we further calculate the Gegenbauler 
moments, and establish a model for kaon twist-3 wave- 
functions with the help of BHL prescription, which 
have a better endpoint behavior and shall be helpful for 
estimating the kaon involved inclusive or exclusive pro- 
cesses. 
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TABLE VI. The normalization parameters and the first two moments of <j>p :tT and their corresponding Borel windows. 



M 2 (GeV 2 ) 
( M ?) 2 (GeV 2 ) 


[0.799, 0.976] 
1.068 - 1.280 


Af 2 (GeV 2 ) 


[0.895, 1.066] 
0.539 ±0.012 


Af 2 (GeV 2 ) 


[0.833, 1.109] 
0.580 ±0.010 


M 2 (GeV 2 ) 
(^)(GeV 2 ) 


[0.974, 1.052] 
1.229 - 1.316 


M 2 (GeV 2 ) 


[0.831, 1.163] 
0.313 ±0.012 


M 2 (GeV 2 ) 


[1.095, 1.340] 
0.248 ± 0.004 



As a final remark, by setting the current quark mass 
m a = 0, we can obtain the sum rules for the pion distri- 
bution amplitudes <pp a , whose normalization parameters 
and the first two non-zero moments together with their 
corresponding Borel windows are presented in Tab lVIl 
By adding all the uncertainties in quadrature, we obtain 
*4|iGeV - 1.104±g;gtg GeV and ^liGcV 



1 14Q+0 033 



GeV. If taking the ct s correction into consideration, which 
increases the leading-order results by 15% — 20% [46|, we 
shall obtain ^liGcV € [1.20,1.38] GeV and nZ\iGeV £ 
[1.28, 1.421 GeV. 



1 



<(o)^£T^(o) 
C(o)^^V^(o) 



+j l x ll x v x p (Q i>%(0)D r D D r ipp(0) Oj . (Al) 

The results for the first and the second terms are well 
known, 



(o|C(o)^(o)|o> = -(^>^ 



and 
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As for the third term (0 



C(o)^V^(o) 



pyvj Oj, basing 

on its color and Dirac-gamma structures, it can be rewrit- 
ten as 



Appendix A: Details for the formulas of 

(0|&(z)$(i/)|0) and (0\fc(xW b p(y)G*,\0) 

Under the background field approach, ip(x) can be ex- 
panded around x = [24| . 



Xjj(x) = ^(0) + V(0)^"x Q + -ijlQffiPff X a Xp 

X^X^X'-y I ' ' ' . 



± 



Then, the matrix clement ^0 
panded around x = as 

(0|C(^(0)|0) 

= (o|c(o)^(o)|o> + ^(o 



r a (x)^(o) 



) can be ex- 



C(o)^(°) 



Utilizing the equation of motion of the background 
quark field and the equation [D** ,D] = -ig s T A G A ^ u , 
" (W) + M (9*#°TGil>) and 



we obtain C 
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D = -^g (g s ipaTG?p) : where {g s il)aTGijj) is the ab- 
breviation of (g s ^a l _ lu T A G A ^ v il)). The fourth term 



o ^(o)^tr£>^(o) 



0) can be treated similarly, 



5 ba (EYg up + Fl v g vp + GhfgT)p a 



with E 

F = 



G 



<)(ix:s 



96x3 
4? m 



^(g^aTG^) and 
gg^g (g s 4>aTGip). Taking use of 



translation invariance of matrix element, we finally ob- 
tain 



(0\r a (x)^(y)\0)=5 ba ^) 

+ (g^aTGiP) 



1 im 
12*** + 48 ( ^ 



1 



96 x 2 



- V) 2 9f3a + 



)*.- gg(*-w) ^"9676 
(x - y) 2 {^- tf)p 



(x - yf{^t- $)fs a 



96 x 9 



(A2) 
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The matrix element ( 



V£(aO^(y)G£, 



be expanded it around x = 0, 



0) can also Using the equation of motion of the background quark 
field together with the following equation [471 ]: 



(0\r a (x)^(y)G^\0) = 

(o \r a (o)4(o)G$ u \ o) + Xp (o |c(o)^ p g^( )| o; 

+2/p (o|c(0)G^^(0)|o). (A3) 
Obviously, 



(0 |C(0)^(0)G^| 0) = — (jaTGyj) {o,„) Pa (T A ) b 



Utilizing the equation [T A G A V , < D P ] = T A G A „. P , one can 
derive 





-(0 



<(0)^G^(0) 
C(0)G^^ P ^(0) 0) - (0 |C(0)^(°) G ^;1 °> 



(0|C(0)^(0)G^/|0> = 



432 ' 



we obtain, 



C(0)^ P G^^(0) 



1 

864 
im 



{e p ^r)0a{T A ) ha 



and 



C(0)G^^(0) 



= 



1 

864 

" n (^TG^) + ^-g s (^) 2 



9s (-0-0) 2 {g P niv - g P u r r P )p a (T A ) ba - 

(e Pflv ^ 5 r)0a(T A ) ba . 



384 X1 r/ 864 
Then, we finally have 



1 



(0 \r a (x)4(y)G A „\ 0) = — (^aTGip) {a^)p a {T A ) ba + \ -—g s ( 9pp 1» - g^M* + v) P 



192 

+i{x-y) p 



^TG^) + ± 9s (W)» 



(e ww 757")/3a (^) 



(A4) 



It is found that Eqs. (|A2|A4|l agree with those of Ref . |48| / 



(Eqs.(22,29) there), except that for (0 ip^(x)ipjj(y) 



) other than (y — x) f 



ip%,(x)ipp(y)G A v 0) the last term should be (a; - y) f 



there is no dimension-six term and the coefficient before 

and for 



, {gsipaTGijj) should be 96 ' x9 other than 



96x12 
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